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Abstract 

We study the spectrum of scalar glueballs in SQCD-like theories 
Oh| whose gravity description is in terms of N c D5 color branes wrapped 

(— i ' on an S 2 inside a CY3-fold, and Nf backreacting D5 flavor branes 

wrapped on a non-compact two-cycle inside the same CY3-fold. We 
show that there exists a consistent truncation of the ten-dimensional 
Type IIB supergravity system to a five-dimensional non-linear sigma 
model consisting of four scalars coupled to gravity. Studying fluctua- 
tions of the scalars as well as the metric around particular backgrounds 
allows us to compute their spectra. A few different backgrounds share 
CN ' the same qualitative features, namely that the mass of the lightest 

scalar glueball increases as the number of flavors is increased, until 
one reaches the point Nf = 2N C after which the opposite behaviour is 
obtained. We show that the five-dimensional non-linear sigma model 
obeys Seiberg duality, and demonstrate this explicitly for the spectra 
rN ! of a class of backgrounds that are Seiberg dual to themselves. 
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1 Introduction 



Since the AdS/CFT conjecture was originally proposed in [lj and sub- 
sequently refined in [2J [3] , it has been extended to apply to more gen- 
eralized settings with less supersymmetry and also backgrounds that 
do not asymptote to AdS in the UV. In this paper, we will study 
the spectrum of scalar glueballs in SQCD-like theories, whose grav- 
ity description is in terms of N c D5 branes wrapping an S 2 inside 
a CY3-fold, and Nf backreacting D5 flavor branes wrapping a non- 
compact two-cycle inside the same CY3-fold. The dual field theory 
is believed to be similar in the IR to N = 1 SQCD with a quartic 
superpotential for the quark superfields [4j. However, the full the- 
ory cannot be dual to SQCD for a number of reasons. It does not 
have an SU{N f ) x SU(N f ) x U(l) R global symmetry as SQCD does, 
but instead only one SU(Nf) (broken further to U(l) N f by smearing 
the flavor branes as will be discussed later). Also, for Nf < N c , the 
Affleck-Dine-Seiberg superpotential [5] tells us that SQCD does not 
have a vacuum, whereas for the systems we will study backgrounds 
exist with Nf < N c . The glueball spectra of the theories that we will 
study in this paper are not known since before. Using holographic 
techniques, we will find how the mass of the lightest scalar glueball 
in the spectrum depends on the number of flavors for a few different 
backgrounds. 

The holographic prescription for computing the glueball spectrum 
is to study fluctuations around a particular background and look for 
solutions that satisfy correct boundary conditions in the IR and UV. 
This is only possible for specific values of K 2 = — M 2 , where K is 
the four-momentum of the fluctuations, and these special K 2 give us 
the spectrum. Technically, it is difficult to carry out such a study in 
ten dimensions. However, we will show that there exists a consistent 
truncation to a five-dimensional non-linear sigma model consisting of 
four scalars coupled to gravity. This will simplify the analysis. In [6], 
an explicitly gauge-invariant formalism was developed for studying 
such systems, and formulas for the linearized equations of motion for 
the fluctuations were given in terms of a superpotential W. Here, we 
will derive the generalized versions of these formulas, which hold for an 
arbitrary potential V not necessarily obtainable from a superpotential. 

In the gravity picture, Seiberg duality is realized for these theories 
as a diffeomorphism, i.e. just a change of variables [U [TJ. There- 
fore, the background itself does not change under Seiberg duality, but 
since we have changed variables, the dictionary interpretation of the 
dual field theory is changed. We show that the Lagrangian of the 
five-dimensional non-linear sigma model is invariant under a set of 
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transformations of the scalar fields and N c —> Nf — N c . It follows 
that anything that can be computed within this framework will obey 
Seiberg duality. 

The backgrounds correponding to the setup described above that 
have been found fall into two categories known as Type A and Type N 
[3 IE]- Type A backgrounds are special cases of Type N backgrounds 
for which the VEV of the gaugino condensate as well as the mesons are 
zero. In this paper, we will study the spectrum of a few backgrounds 
of Type A for which the dilaton grows linearly in the UV. In the IR, 
there are different possible behaviours for the background (known as 
Type I, II and III [8]) corresponding to different vacua in the dual 
field theory. These backgrounds have a singularity in the IR which 
is "good" according to the criterion given in [10J, and are believed to 
capture the non-perturbative physics of the dual field theory. This cri- 
terion states that the goo component of the metric should not increase 
as we approach the singularity (the idea is that proper energy exci- 
tations should correspond to lower and lower energy excitations from 
the point of view of the field theory as one approaches the singularity 
in the IR). 

The D5 flavor branes are smeared along the transverse angular 
coordinates, breaking the SU(Nf) global symmetry to U(l) N f (this 
procedure was first introduced in the context of flavor branes in [9]). 
The consistent truncation to five dimensions does not contain fluctu- 
ations of the gauge fields on the branes. However, it still contains 
fluctuations of the Ramond-Ramond 3-form Fr$\. Therefore, when 
Nf ~ N c , the fluctuations that we consider mix glueballs and mesons. 
Since the fluctuations do not involve the gauge fields on the brane, 
the meson-glueballs whose spectrum we compute are [/(l)^-singlets. 

Imposing the boundary condition on the fluctuations in the IR 
that their kinetic terms are regular, and in the UV that the fluctu- 
ations correspond to normalizable modes, we find that the mass of 
the lightest scalar glueball increases as the number of flavors is in- 
creased, until the point Nf = 2N C is reached after which the opposite 
behaviour is observed. For a particular class of backgrounds that are 
Seiberg dual to themselves, we demonstrate explicitly that the spec- 
trum obeys Seiberg duality. 

There is by now a large literature on systems with back-reacting 
flavors. In the future, it would be interesting to apply the same tech- 
niques to study the glueball spectra of the various systems studied in 

[Zl IB1 [ID H21 U31 US HS1 031 UB1 US1 12D1 12U [221 [221 123 12SI I2ZI 
12312313211311132113311311133133133 ESI- 

This paper is organized as follows. In section 2, we describe the 
general setup and the backgrounds that we will study. In section 3, 
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we introduce the five-dimensional non-linear sigma model, discuss the 
Seiberg duality it obeys, and present the holographic techniques that 
we will use to compute the spectra. Section 4 contains the compu- 
tation of the spectra. Finally, we summarize our results in section 
5. 



2 Gravity Duals of SQCD-like Theo- 
ries 

The backgrounds we will be interested in are obtained from wrapping 
N c D5 color branes on an S 2 inside a CY3-fold, then adding Nf flavor 
branes that wrap a non-compact two-cycle inside the same CY3-fold. 
This is described in detail in |4], where evidence is given for that the 
backgrounds obtained are dual to a field theory with similar behaviour 
in the IR as M = 1 SQCD with a quartic superpotential for the quark 
superfields. 



2.1 Action and Equations of Motion 



We will now write the Type IIB supergravity action and the equations 
of motion that follow from it. The action (in Einstein frame) is given 

by 



S = S IIB + S iflavors \ 



(1) 



where Sub describes Type IIB supergravity in the truncation to the 
metric, the dilaton, and the RR 3- form (g^ u ,(p,F^), and 5 , (/ /avors ) jg 
the action of the flavor branes. We have that 



Sub 



1 



2k- 2 



d w x, 



-(J 



R 



-e*FL 
12 ( 3 ) 



(2) 



Choosing coordinates p, 6, ip, 6, tp, ip), the flavor branes extend along 
the external coordinates x^, the radial coordinate p, and the angular 
coordinate ip. Their action is given by 



s (flavors) _ Tm ^ 



[ d 6 xe^ 2 v /=^+ / P[C 6 ) 
JMe J Me 



(3) 



where <7(g) is the determinant of the pullback of the metric to A4q, the 
world volume of the flavor brane, and similarly P[Ck] is the pullback 
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of Cq. After smearing along the angular coordinates (9,<p,6,<p), this 
becomes 



g(flavors) 



T D5 N f 
(4tt) 2 



where 

S7 4 = sin sin 9d9 Ad9 A dip A d^. 
The equation of motion for the dilaton is 

=0/2. 



(4) 



(5) 



1 



1 



-5 



\/~3(io) 



sin sin # = 0, 



(6) 



while Einstein's equations read 



1 



1 



+ 



12 



(7) 



where 



N f 



^ sin sin fle*/ 2 »S "^=^ ■ 

! - ,/ "'"" r! <S / 9y ( 6 ) y/-p(l ) 



(8) 



Finally, Maxwell's equation for the is given by 

(V=ge*F^ x ) = 0. (9) 



2.2 Type A Backgrounds 

In this paper, we will be interested in so-called Type A backgrounds. 
For these backgrounds the VEV of the gaugino condensate is zero, 
as are the VEVs of the meson matrix. Type A backgrounds can be 
obtained from starting with the ansatz [3] 



ds 2 =/ u 2 e 2 ^ 



V~ 2 dxi t3 + e 2k dp z + e 2h (d9 2 + sin^ 0ckff)+ 

2k 



(d9' z + sin 2 9d(p 2 ) + ^—(dij) + cos 9 dip + cos 9d<pf 



(10) 



6 



r ^ sin QdO A dip + — - ^ c sin 8d9 A dtp 
4 4 
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(dtp + cos 6*(i</3 + cos 9 dip) , 

with /x 2 = a'g s . Here, /, k, h, and g are taken to be functions of the 
radial coordinate p. For the backgrounds that we will be interested in 
the IR is at p = and the UV at p = oo. 

Making a change of variables to P, Q, and Y through 

r 2h_ P + Q 

4 ' 

e 2 ~9 =P-Q ( 12 ) 
e 2k =4Y, 

the BPS equations can be solved as (/ = 0/4) [8] 



Q =Q + (2iV c - N f )p, 



(13) 



(P 2 -Q 2 )Y : 



where cfto and Qq are integration constants, and P satisfies a second 
order differential equation given by 

» / / s fP' + Q' + 2N f P'-Q' + 2N f \ 

2.3 IR and UV Expansions 

We will be interested in backgrounds for which P grows linearly in the 
UV0 For Nf < 2N C , these have the UV expansion (around p = oo) 
given by 

P UV = {2N C - N f )p+ (1 + Q )N C + 4(2 ^ C jV/) P' 1 + O (p~ 2 ) , (15) 



whereas if Nf > 2N C 



Nf(N f -N e ) „_ ll/n ,„_ 2 , 



P UV = -(2N c -N f )p-((l+Qo)N e -N f )- 4( J 2 ^_ jfjf P~ l +° (P 



(16) 



1 The alternative is that the P grows exponentially in the UV, in which case the spec- 
trum only contains a continuum. 
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1 2 3 4 5 P 

Figure 1: P as a function of p for the Type A background with Type II IR 
behaviour and Qq = 1.2. The different lines correspond to different number 
of flavors: dotted is N f = N c , dashed is N f = 1.4JV C , and solid is N f = 1.8N C . 




Figure 2: P as a function of p for the Type A background with Type III 
IR behaviour. The different lines correspond to different number of flavors: 
dotted is Nf = N c , dashed is N f = IAN C , and solid is N f = 1.8N C . 




1 2 3 4 5 P 

Figure 3: P as a function of p for the Type A background with Type II IR 
behaviour and Q = 1.2. The different lines correspond to different number of 
flavors: dotted is N f = 2.2N C , dashed is N f = 2.QN C , and solid is N f = 3N C . 
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Figure 4: P as a function of p for the Type A background with Type III 
IR behaviour. The different lines correspond to different number of flavors: 
dotted is N f = 2.2N C , dashed is N f = 2.6N C , and solid is N f = 3N C . 

In the IR, there are several different possible behaviours for the 
background. In this paper, we will focus on two different ones, that 
of Type II and Type III [8j. For Type II and Qq > we have that @] 

pM = Q + 4h 1 ^- (w e + || + JV>) p + O (> 2 ) . (17) 

There are two integration constants: Qq and h±. Solutions exist that 
interpolate smoothly between the Type II IR and the linear dilaton 
behaviour in the UV [3]. In order to obtain such solutions, one must 
dial the integration constants Qq and h\, essentially making hi a func- 
tion of Qo- This leaves us with one free parameter Qq for the Type II 
solutions. 

In the case of Type III, Qq = and P has the following behaviour 
in the IR: 

P W = 4 V V 3 - + ^V /3 + O ( P ^) . (18) 

Now, requiring that the solution has the UV asymptotics of Pjjv com- 
pletely fixes the one parameter h±. Both the Type II and Type III 
backgrounds have a singularity in the IR, which satisfies the criterion 
for being a "good" singularity given in |10| . 

Figure [1] shows P as a function of p for the Type A background 
with Type II IR behaviour and Qo = 1.2 for a few different number 
of flavors Nf < 2N C . Figure [2] shows the corresponding plots for the 
Type A background with Type III IR behaviour. Figures [3] and H] are 
the same as Figures [T] and [5] but for flavors Nf > 2N C . 
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2.4 Dual Field Theory 

We will now describe in more detail some aspects of the field the- 
ory conjectured to be dual to the backgrounds studied above. First, 
consider the case of no flavors. Then the dual field theory is a four- 
dimensional M = 1 super symmetric field theory, obtained from a 
twisted compactification of six-dimensional SYM on S 2 where the 
twisting is such that it preserves four supercharges HQ] . At weak 
coupling, this field theory consists of a massless vector multiplet, V, as 
well as a Kaluza-Klein tower of massive chiral and vector multiplets, 
<& k and V k . The infinite number of KK modes reflect the fact that the 
UV completion is not given by a quantum field theory (in fact, it is 
given by a Little String Theory). If it were possible to separate the 
scale set by the size of the S 2 from the scale A at which the theory 
becomes strongly coupled, we would have a gravity dual of M = 1 
SYM. Unfortunately, this is not the case. The Lagrangian of the field 
theory without flavors has the generic form 

C= I d^Y.^y^k + m k \V k \ 2 ) + 
k 

W a W a + Y (Wk, a W? + n k |<D fc | 2 + W(<f fc , V k )) , 

k 

(19) 

where W a and W kja are the curvatures of V and V k , and m k and /i k 
are the masses of the massive vector and chiral multiplets comprising 
the Kaluza-Klein tower. The superpotential W(& k , V k ), governing the 
interactions between the KK chiral and vector multiplets, is given by 

W(* k , V k ) = Y + Y f(*k)W kia W£. (2Q) 

i,j,k k 



With the introduction of flavors, we also have to add the terms [1] 




to the Lagrangian. Here, a, b = 1,...,N C are indices of the funda- 
mental and anti-fundamental representations of SU(N C ) and i,j = 
1, . . . , Nf are indices of the fundamental and anti- fundamental repre- 
sentations of SU(Nf). Since the smearing procedure described above 
breaks SU(Nf) to U(l) N f , the Kp must serve the role of breaking this 
symmetry in the field theory, however its exact form is not known. In 
principle, we could have also considered the more general case where 
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there is a superpotential for the flavors too. Again, the form of such 
a superpotential is not known. 

As mentioned above, it is not possible to separate the scale set by 
the size of the S 2 from the scale A at which the field theory becomes 
strongly coupled. Nevertheless, we can imagine integrating out at 
least some of the KK modes, if not all the way down to A. This then 
gives rise to an effective superpotential W e ff containing quartic terms 

Wef f ~J2%r(Q*Qj) 2 - (22) 

P,h3 P 

This means that in the IR the field theory is similar to M = 1 SQCD 
with a quartic superpotential (bearing in mind that not all the KK 
modes can be integrated out in this fashion). 



3 5d Formalism 



3.1 Non-linear Sigma Model 

There exists a consistent truncation of the lOd system to a 5d non- 
linear sigma model consisting of four scalars coupled to gravity. The 
action of this 5d model is obtained by plugging in the ansatz ()10p into 
the Type IIB supergravity action and integrating over the angular 
coordinates (see Appendix A for a detailed derivation). After making 
a change of variables as 



(23) 



A = ^(8f + 2g + 2h + k) - log 2, g = -g + h + log 2, 

p = - \ (4/ + ~g + h + 2k) + \ log 2, x = 2/ + g + h - log 2, 
b 2 

and also changing to the radial coordinate as dz = e A+k dp, one obtains 



C 



e// _ 4^(47r) 



R 



\G ab d<$> a d<s> b - 



-5d - — ~ v ->/ 

where the non-linear sigma model metric G a b is given by 
Go 



<!ia 8 ( 2' 6 ' 1 ' 4 



(24) 



(25) 



and the potential V is given by 



V 



-2(g+2(p+x)) 



128 



16 (_4e 9+6p+2x (1 + e 29 ) + e 4s + l) + 



(26) 
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The backgrounds we will be interested in are of the form 



ds 2 = dz 2 + e 2A dx{^ (27) 

where A is a warp factor. 

We note that V can be written in terms of a superpotential W as 
follows: 

V = \w a W a - ^W 2 , (28) 



with 



(29) 

4(1 + e 29 + 2e 9+6p+2:E ) 



where W a = dW/d$ a , and indices are lowered and raised using the 
non-linear sigma model metric G ab and its inverse. The first order 
equations of motion following from the superpotential are 

A' = -™ 

3 ' (30) 

$ /a = W a , 

and agree with the BPS equations given in Appendix B of In 
Appendix D, we show how to derive the expression for W from the 
BPS equations. 

The equations of motion for the scalars following from the 5d action 
read 6j 

V 2 $ a + g a bc g^(d^ b ){d u ^ c ) - V a = 0, (31) 
whereas Einstein's field equations read 

-iV + 2G ab (d^ a )(d^ b ) + \g^ V = 0- ( 32 ) 

Here the tilde is used to denote a five-dimensional quantity, and Q\ c 
is the Christoffel symbol with respect to the non-linear sigma model 
metric. One can verify that every solution to these 5d equations of 
motion also solves the full lOd Type IIB supergravity equations of 
motion for the particular choice of non-linear sigma model metric G ab 
and potential V above. This shows that the five-dimensional non- 
linear sigma model is a consistent truncation of the lOd Type IIB 
supergravity system. 
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3.2 Seiberg Duality 

For the models considered in this paper, Seiberg duality is realized on 
the gravity side as a diffeomorphism. While the background does not 
change, the change of variables means that the dictionary describing 
quantities in the dual quantum field theory changes. 

In terms of the variables (P, Q, Y, (ft), Seiberg duality transforms 

in 

W (33) 
N c — > Nf — N c , 

leaving P, Y, and (ft unchanged. Using the relations 

e^(P 2 -Q 2 )VY 



e 3A 



e 6p 



e 2x 



16 

P + Q 
P-Q 1 

4 e -<t> 

\yp2_Q2y 1 

e*{P 2 -Q 2 ) 
16 



(34) 



we see that in terms of the 5d variables, a Seiberg duality simply 
takes the form g — )■ —g. It is straightforward to see that both the 
potential V and the non-linear sigma model metric Gij are invariant 
under this transformation. It follows that the whole 5d theory exhibits 
Seiberg duality, and therefore anything that we can compute within 
this framework, including the spectrum, will manifest Seiberg duality. 

Considering that Seiberg duality is normally only a duality in the 
IR, it may seem odd that the whole 5d Lagrangian is invariant under 
the Seiberg duality transformations. However, in [51] it is argued that 
M = 1 SQCD with a quartic superpotential for the quark superfields 
can satisfy an exact Seiberg duality, where not only the IR of two differ- 
ent field theories are the same, but they in fact represent two different 
descriptions of the same renormalization group flow. The fact that 
in our setup Seiberg duality corresponds to diffeomorphisms supports 
the view that Seiberg duality is exact for these backgrounds. However, 
it is not clear how this happens from a field theory perspective. One 
can make the following schematic argument in the IR. Starting with 
a superpotential W = /i/x -1 (QQ) 2 , we first Seiberg dualize, obtaining 
W = h/j,M 2 + qMq. For Nf < 2N C , h is a relevant operator, so we can 
integrate out M, solving = dW/dM = 2hfLM+qq. Plugging M back 
into W gives us an effective superpotential W = —h~ l n~ l (qq) 2 /2. As 
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can be seen, this superpotential is of the same form as the one that we 
started with, but describing a theory with N c = Nf — N c colors and 
Nf flavors. Furthermore, the coupling h has been inverted, consistent 
with that for Nf > 2N C the coupling appearing in front of the quartic 
term of the superpotential is irrelevant. It is, however, clear that this 
argument only works in the IR. In order to argue from the field the- 
ory that Seiberg duality is exact, one would have to take into account 
the KK modes that become important in the UV. From this point 
of view, it still remains somewhat mysterious why the backgrounds 
we are considering seem to have an exact Seiberg duality and what it 
means in terms of the dual field theory. 

3.3 Holographic Techniques for Computing 
Spectra 

The holographic prescription for computing the spectrum of glueballs 
is to study flucuations around the background, expanding the equa- 
tions of motion to linear order and look for solutions for the fluctua- 
tions that satisfy the correct boundary conditions in the IR and the 
UV. This is only possible for special values of K 2 , where K is the 
space-time momentum of a plane wave. These K 2 correspond to poles 
of the correlator (OO) (where O is the operator in the dual field theory 
corresponding to the fluctuation in question), and give us the scalar 
glueball spectrum of the dual field theory. In order to find them, we 
will employ a generalization of the numerical method described in [12] 
that in effect evolves solutions from both the IR and the UV, then 
determines whether they can be matched up smoothly at a midpoint 
(see Appendix E for a detailed description of the numerical methods 
used) . 

In [B], an explicitly gauge invariant formalism was developed to 
compute spectra holographically. Formulas for the linearized equa- 
tions of motions for the fluctuations were given in terms of a superpo- 
tential W. Here, we will give the generalized corresponding formulas 
that are valid for any potential V. Following J6], we expand the metric 
and the scalar fields in fluctuations that are gauge invariant. First, we 
write the 5d metric in terms of a lapse function and shift vector n l as 

ds 2 = {rtrii + n 2 )dz 2 + 2n i dx t dz + g ij dx l dx j . (35) 

Then, we expand to linear order in the fluctuations (o a , b, c, J , tij) 
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using the prescription given in [6]: 

_> $ a + a a , 
n-> 1 + b, 

i !U/ i # . (36) 
n l -> e 2A (^ + -c), 

Oil , V 

.'/-;, ' ' ('/>./' • c,.,)- 

Note that this is not just a probe calculation (in which the fluctuations 
of the metric are put to zero), but both the scalar fields and the metric 
are dynamical. 

For a general potential V, the linearized scalar equations of motion 

are 



Dl + 4A'D, + e' 2A U 



a a -(V a lc -K a bcd <S>' b <t>' d K- 
& a (c + d z b) -2V a b =0. 



Here, V a = dV/d$> a , and, as before, indices are lowered and raised 
by the sigma model metric G ab and its inverse. If an index is placed 
after a "|", it means that the covariant derivative with respect to G ab 
should be taken: 

A a \ b = D b A a = d b A a - G c ab A c , (38) 

where 

<5\c = \G ad {d c G db + d b G dc - d d G bc ). (39) 

The Riemann tensor with respect to the non-linear sigma model metric 
is given by 

= 9 c G a bd - d d Q\ c + Q\ e Q\ d - Q a de Q\ c . (40) 
Also, a "background-covariant" derivative D z has been defined through 
D z <p a = d z ip a + Q\^'\ c . (41) 
Similarly, from the linearized Einstein's equations, we obtain 

6A'c + 4*aP*0 - W a a a - 8Vb = 0, 

1 (42) 
--no i + 3A'd i b-2& a d i a a =0, 

where prime denotes the derivative with respect to the radial coordi- 
nate z, and &' a = G ab & b . 
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As in [6], (|42p can be solved algebraically in momentum space 
giving the metric fluctuations b, c, and in terms of the scalar fluc- 
tuations a a . Plugging back into (f3T|) then gives us a linear differential 
equation for just the scalar s: 



D 2 + AA'D Z - e~ 2A K 2 



,ra */6*/«i 4($ ,a U c + V a ® c ) 16U$' a $ 



(43) 



0. 



3A' ' 9A> 2 

In the special case where V can be written in terms of a superpotential 
W, this formula agrees with the one given in [6]: 

(44) 

a c = 0. 



btD z + Wi 



\b 



Changing the radial coordinate as dz = e A+k dp, ()43[) becomes 



5 b a d 2 + S%d p + T£ - 5 a b e 2k K 2 



0, 



(45) 



with 

S a b =2Q\ c d p & + 4 (0 pP + d p A) 51 

T b a =d b g a cd d p ^ c d p ¥- 

' 4(V a d p $ c + V c d p <S> a ) 



4e -8 P 



3d p A 



+ 



16Vd p $ a dp§> c 
9(d p Af 



G rh + d h V a 



(46) 

This is the second order linear differential equation for the scalar fluc- 
tuations that we need to solve for different values of K 2 imposing the 
correct boundary behaviour in the IR and UV. In the IR, we will re- 
quire that the kinetic terms for the fluctuations are regular. In the 
UV, we will require that the fluctuations are normalizable. 



4 Scalar Spectra 

In this section, we will study some different Type A backgrounds, 
and compute the mass of the lightest scalar glueball as a function of 
the number of flavors. First, we will work out what the boundary 
conditions for the fluctuations are in the IR and UV. Then we will 
solve (|45p numerically for different values of K 2 using the methods 
outlined in Appendix E. This gives us the spectrum. In the following, 
we will put Xf = Nf/N c , and rescale P — > N C P and similarly for Q 
and Y. All masses given are in units of a'g s N c . 
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4.1 Boundary Conditions in the UV 
4.1.1 N f < 2N C 

We will now expand the differential equations for the scalars (|45[) in 



the UV. For Nf < 2N C , the background is given by (fT5|) . We obtain 
that 



S a b =25 a b +0{p- 1 

( -4 I 

T = 



\ 



6 -1 
6 -3 
12 2 



"2 \ 

l 

2 

-3 y 



-IN 



(47) 



A basis that diagonalizes T to leading order is given by 

\ 



B 



( "I 

l 

2 
1 

V o 





1 



1 



(48) 



•1 / 



such that B~ 1 TB is diagonal. To leading order, this diagonalizes the 
differential equations for the fluctuations, and we obtain four indepen- 
dent differential equations 



d 2 a 1 + 2d p a 1 
dW + 2d a 2 



+ K 2 )a 1 



p 1 

d 2 p a 3 + 2d p a 3 



d 2 a 4 + 2d D a 4 



K 2 )a 2 
-K 2 a 3 
-K 2 a A 



--{) 
--{) 
--{) 
--{) 



(49) 



with solutions a 1 ' 2 ~ 6 (-i±V9+^)p and a 3,4 _ e (-i±vTW)p. Note 
that if we imagine expanding the fluctuations as a a = \ a {p)e^~ l ^ Ca ~ M1 ' > 
with X a (p) = Yl n a nP ba ' n i the above analysis captures the C a but not 
the function X a (p). Therefore, the exponential factors are in general 
multiplied by powers of pE 

However, the exponential behaviour is all 
we need for setting up the boundary conditions in the numerics. We 



2 The validity of the expansion in powers of p^ 1 hinges on that it is possible to find a 
basis in which the components with different exponential behaviour do not mix. In the 
case of [45], this can be checked explicitly. However, in that case, P is exponentially close 
to P = 2N c p in the UV, so that we can always work with analytical expressions. In the 
present case, only the UV expansion of P in powers of p~ x is known. While we cannot 
verify that there exists a basis in which the different exponential behaviours do not mix, 
the fact that we can go to reasonably high cut-offs in the UV (around p = 15) before the 
numerics break down suggests that such a basis exists. 
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are interested in the subleading behaviour so we pick the minus signs. 
In [12], a normalizability condition for the fluctuations was given: 



dze 2A G ab ^ a ^ b = I dpe M+K G ab ^ a ^° < oo. 



3A+kri „/,a„/,& 



(50) 



In our case, we have in the UV that (xf = Nf/N c ) 



e 3A+k = e 2p+2<j, 



(51) 



so that the subdominant fluctuations are always normalizable, while 
the dominant ones are not. Let us also point out that for M 2 > 1 or 
M 2 > 9, we start getting oscillatory behaviour for fluctuations in the 
UV, signalling the start of a continuum. 



4.1.2 N f > 2N C 

Similar considerations as in the last section (but with a different B) 
now le ad to solutio ns of the form a 1 ' 2 ~ e (~ 1± V^ : i^f-^)^)p anc } a 3 > 4 ^ 
e (- 1 ±\/ 1 +( x f- 1 ) K2 )p (times powers of p). Note that the appearance of 
xj is Seiberg duality at work (restoring units, it takes g s a'N c K 2 — > 
g s a'N c (x f -l)K 2 ). 



4.2 Boundary Conditions for Type II in the 
IR 

For Type II backgrounds, it is natural to expand the fluctuations in 
the IR as 

oo 

a a = J2<P n/2 - (52) 

n=0 

We will choose boundary conditions such that the kinetic terms of the 
scalars do not blow up, i.e. that the derivative d p a a does not blow 
up in the IR. This fixes af = and, after plugging in the ansatz (|52[) 
into the differential equations for the fluctuations (|45p . leads to four 
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linearly independent solutions 



0(1) 



0(2) 



0(3) 



0(4) 










+ 





w 





1 





W 



+ 



(-2 + 



2xj_ Whj \ 

Qo ^ ~Ql ' 

2+2Q -x f 

2Qo 
2+AQo-x f 

2Qo 
2-Xf 



\ Qo 

Qo 
at 5x f 

(i • 2(77 „, 
+ ^ + A 



/ 

Qo 
3 



+ 7? 



2Qo 
3a; / I 6 



P + 0(P 3/2 ), 



(53) 












1 


+ 


w 




/ON 












+ 



Qo " 
x f 
4Qo 



Qo 



+ 



2hf 



-A _|_ ±1_ _|_ 
^ 4Qo ^ 2A 

i 1 
2Q ^ 77? 



0(P 3/2 ), 



i. \ 

l 



2Qo 
U_ , 
8Qo ^ 

8Qo ^ 4/i? 
"1L + l 



4hf 



p + 0(p 3 / 2 ). 



\4Q ^ 277?/ 

This fixes our boundary conditions in the IR. 



4.3 Boundary Conditions for Type III in the 
IR 

For Type III backgrounds it is natural to expand the fluctuations as 

oo 

a a = Y,<P n/3 - (54) 

n=0 

We see that the requirement that the derivatives of the fluctuations 
do not blow up in the IR now leads to of = = 0, which is a stronger 
requirement than for Type II, and consequently leads to fewer than 
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Figure 5: The mass squared of the lighest scalar glueball as a function of 
the number of flavors for a couple of Type A backgrounds with Type II IR 
behaviour: Q = 20 (squares) and Q = 1.2 (dots). 



four allowed linearly independent solutions in the IR: 



a (i) 



1(2) 
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V i , 




W 






2 


) 



P i/3 + 0(p 



5/3 \ 



n 4 /3 



(55) 



Now that we have fixed the boundary conditions, singling out a 
number of allowed linearly independent solutions in the IR and UV, 
the question becomes whether for a particular value of K 2 it is possible 
to find linear combinations of the allowed solutions in the IR which 
when evolved towards the UV can be written as linear combinations 
of the allowed solutions in the UV. The numerical methods used for 
determining this are outlined in Appendix E. 



4.4 Results 

Figure [5] shows the mass squared of the lightest scalar glueball as a 
function of Xf = Nf /N c for a couple of Type A backgrounds with Type 
II IR behaviour. As can be seen, the mass increases with the number of 
flavors, until the special point Nf = 2N C (where the theory has certain 
peculiar properties) where it reaches the start of the continuum, M 2 = 
1, and after that decreases as a function of the number of flavors. 

Figure [6] shows the mass squared of the lighest scalar glueball as a 
function of xj = Nf/N c for the Type A background with Type III IR 
behaviour. Again the same pattern can be seen. 
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Figure 6: The mass squared of the lighest scalar glueball as a function of the 
number of flavors for the Type A background with Type III IR behaviour. 



■ 



Figure 7: The mass squared of the lighest scalar glueball as a function of the 
number of flavors for the Type A background with Type III IR behaviour 
with the Seiberg dualized spectrum superimposed. 



Under Seiberg duality, the integration constant Qq — > —Qq. Since 
for a Type A background with Type III IR behaviour Qq = 0, such 
backgrounds are Seiberg dual to themselves (Q = (2N C — Nf)p — > 
(2N C — Nj)p). In Figure the dots are the same as in Figure El 
and the squares are what is obtained under Seiberg duality, mapping 
points as Xf — > f^ . Also, under Seiberg duality, M 2 — > (x/ — 1)M 2 , 
since we are working in units a'g s N c — > a'g s (Nf — N c ). As can be 
seen, the Seiberg dualized spectrum falls on the same trajectory. 



5 Conclusions 

We have been able to find a consistent truncation of the ten-dimensional 
Type IIB supergravity system describing iV c D5 color branes and Nf 
backreacting D5 flavor branes to five dimensions. The five-dimensional 
system is a non-linear sigma model coupled to gravity. In this model, 
Seiberg duality is realized at the level of the Lagrangian, i.e. any 
quantity that we can compute will automatically obey Seiberg dual- 
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ity. 

We have computed the mass squared of the lighest scalar glueball 
for a few different Type A backgrounds, and found that the mass 
increases with the number of flavors for Nf < 2N C , but shows the 
opposite behaviour for Nf > 2N C . For a class of backgrounds that 
are Seiberg dual to themselves, we have seen explicitly how Seiberg 
duality is realized for the spectrum. 

In the future, it would be interesting to apply the same techniques 
in order to compute the spectra of different systems with back-reacting 
flavors. For example, gravity duals that exhibit walking behaviour 
were found in |43|. |4"1"1 14"5] , and in particular one could imagine adding 
flavors to the walking backgrounds of [35] (for which P grows linearly 
in the UV) and study how the spectrum is affected. It would be 
interesting to know what the effect of flavors is on the light scalar 
present for these backgrounds. We leave these questions for a future 
study. 
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A Deriving the 5d Effective Action 

Here, we will derive the 5d effective action of the non-linear sigma 
model for the more general case of Type N backgrounds, of which 
Type A is a special case. We start with the ansatz 



ds 2 =fi 2 e 2 f 



H~ 2 dxl 3 + e 2m dp 2 + e 2h (d8 2 + sin 2 9dip 2 )+ 



e 2g e 2k 

((a)i + ad6) 2 + (<I>2 — a sin 9dp) 2 ) -\ — ^-(^3 + cos 9dip) 2 



fxN c 



{&i + bd9) A (C02 — b sin Odip) A (£3 + cos 9dip)+ 



dy^ A (d^b(-d6 Awi + sin 9 dip A tD 2 )) + (1 - 6 2 ) sin 9d9 Ad(pAu 3 
sin 9d9 A dip A (dip + cos 9d<p) , 

(56) 
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where p? = a'g s , and 

u>i = cos ipdO + sin ip sin 9 dip, 
0J2 = — sin ipdO + cos ip sin 9d(p, 
0)3 =cfa/> + cos 9 dip. 



(57) 



First, let us consider the case when a, b, f, 5, /i, k, and m only depend 
on p. Plugging (|56|) into the Type IIB action given by [1] 



Siod 



2k (io) 



d 10 x^g- 



R--d„ 
2 M 



12 ( 3 > 



(58a) 



T D5 N f 
(4vr) 2 



d 10 x sin 9 sin tf^/ge; - / C 6 A f2 4 ) , 



4 = sin 9 sin 6>e$ A d9 A dtp A dtp, (58b) 
and performing the integration over the angular coordinates yields 
V(4vr) 3 



Swd = \ ^ [ d^x I dpe s f+^+ 2h+k+m (T — V) , (59) 
2k (io) J J 



where 



-2m 



e 2g-2h 
128 



.12 



N 2 e -Af-2h+<t>-2~g 



128 



-6 2 + — 1 

4 16 



$ 2 g' 2 fk' h'k' 1 

— + — + f ti + 1 h f q + -rh q + -k'q 

64 16 J ^2 8 ,/y 4 y 8 y 



(60) 



and 



V 



-2(2(f+h)+k+2~g) 



256 

8e 2(2/+h+2fc+9) a 2 + 8e 4/+2h+6s a 2 + 16 g4(/+A+*) + 
( a 2 _ X ) 2 e 4(/+fc+9) _ 64e 2(2(/+ft)+fc+g)_ 

16 (a 2 + 1) e 2 ( 2 /+^+ 2 s) + i 6e M N 2 + 8(a _ b) 2^H+^. 



e m~g (jvy _ ( a 2 _ 26a + X ) ^2 + 8e |(4(/+ft+*)+*+^) iV - / 



(61) 



Notice that the Wess-Zumino term, whose only effect is to change the 
Bianchi identity of F3, does not appear in this action, from which the 
Einstein, dilaton and Maxwell equations are derived. 
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Let us change coordinates to 

f = A+p-^, 3 = -A- § + log2-p + x, 

2 2 (62) 

h = -A + | - p + x, k = -A + log 2 - 4p, 

with inverse 

A = \{%f + 2~g + 2h + k) - log 2, 5 = -5 + /i + log 2, 

3 ! ! (63) 

p = - - (4/ + g + h + 2k) + - log 2, x = 2/ + 5 + h - log 2, 

and also change the radial coordinate as dz = e A+m dp. This leads to 
S wd = V ^ (47r)3 / d 4 x [ dze 4A (T — V) , (64) 



2k 2 



with 



1 at-2 <j>- 2x 12 12 ±12 

T = 3A' 2 - \e~ 2 ^ 2 - ^ b' 2 -^--3p' 2 -^--^, (65) 
and 

-2(g+2(p+x)) 

v=- 



128 

e 12 P+ 2 X+4 , ^ 2e 2 9{a _ b) 2 + e 4g + ( fl 2 _ ^ + ^ ^ 

2 (a 2 - 26a + l) e^^NfNc + J 2 ** 2 ***!*] + 8e 2 ° +6p+x+ l N f + 

16 (a 4 + 2 (V - e 6p+2x ) 2 - l) a 2 + e 4 ^ - 4 e s+ 6 P+ 2 * (l + e 2 °) + l) 

(66) 

Recognizing that for a metric given by ds 2 = dz 2 + e 2A dx 2 3 the Ricci 
scalar is (up to partial integrations) equal to R = —12 A' 2 , we can 
write this as the action of a 5d non-linear sigma model 



z/t (io) 



4/i 4 (4vr)^ 

- I a x J azy-g — - -^ju^ u- -v - v ) 

(67) 



where = [g,p, x, (/>, a, b], and the non-linear sigma model metric is 
diagonal with entries G gg = \, G pp = 6, G xx = 1, = \, G aa = 

and Gbb 



e -2g jV 2 e -2x+0 

and Gfefe = 
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B Seiberg Duality for Type N 



We can generalize the arguments of section 3.4 to Type N. For these 
backgrounds, it is convenient to go to the variables P, Q, Y, r, and a 
defined through [7] 



e 2h 



If P 2 -Q 2 



4 \ P cosh t - Q 
e 2 s =PcoshT-Q, 



e 2k =4y, 



Psinh 



(68) 



P cosh r — Q ' 



Seiberg duality corresponds to the transformation 

Q -> -Q, 

a -> -a, (69) 
N e ^N f - N c , 

leaving P, Y, (j), and r unchanged. Using the relations 

3A e 2 ^(P 2 -Q 2 )VY 
G 16 



e 



2 q_ P 2 ~Q 2 



(cosh tP-Q) 2 ' 

, a. (70) 

e 2:c 



e^(P 2 - Q 2 ) 



16 

we see that in terms of the 5d variables, a Seiberg duality takes the 
form 



e 9 



a 



b 



e 2 9 + a 2 ' 
a 

e 2g _)_ a 2 ' 



N c — Nf 
N c ^N f - N c . 

Again, it is straightforward to see that both the non-linear sigma 
model metric G a b and the potential V of the previous section are 
invariant under these transformations. It follows that the whole 5d 
theory obeys Seiberg duality. 
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C Linearized Equations of Motion 



Here we will generalize the results of [6] to cases where the potential 
V can not necessarily be written in terms of a superpotential. The 
discussion closely follows that of [6j. 

C.l ADM formalism 

We start by writing the metric on the form 

where the indices i and j run over the four-dimensional space-time, 
and there is also the radial coordinate (index 0). The tilde is used to 
refer to five-dimensional quantities. The inverse metric is given by 



n 2 ( -n l n 2 g^ + n V ) ' ( 



The second fundamental form is 

JCij = nf°- = --^(do9ij ~ Vjnj - Vj-rtj), (74) 
and one can derive the following relations 

pfe _pfe _ n y 

1 r? J 
T%=-d i n + —IC ij , 
n n 

n \ n z J 

f{] =— (don + n^djn + n l n J /Cjj), 
fg =d n k + rfVin* - nV fc n - InK^rt - n k t° 00 . 

C.2 Expanding the Scalar EOMs 

The equation of motion for the scalars following from the 5d action is 
V 2 $ Q + g%^ v (d^ b )(d v $ c ) - V a = 0. (76) 
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Using the relations of the previous section to rewrite this in terms of 
4d quantities, we obtain 

jd 2 - 2n%d t + n 2 V 2 + rtniVidj - (nJCj + d z Inn - n% \nn)d z + 

\nV l n — d z n l + n^Vj-n 1 + n*(n/Cj + <9 2 Inn — n J <9j Inn) <9j j$ a + 



6 i 



be 



(d z ^ b )(d z ^ c )-2n l {d^ b ){d z ^ 



(n 2 g ij + n l n 3 ){d i $ b )(d j <S> c ) 



2 /-tab 



dV 



= 0. 

(77) 



We will now expand around a background (assumed to be dependent 
only on z) to linear order as follows: 



f + a « 5 

n — > 1 + b, 

n* e 2A (^ + £-c), 



3* 



(78) 



□ 



It will be useful that 



nic) -> -a 2 A + - ( + a,* 1 + - d z z) 



and 



nIC] 



-dd z A + c, 



(79) 



(80) 



where d + 1 is the number of space-time dimensions (so that in our 
case d = 4). At zeroth order, we have 



-dA'<Z>' a - G a bc 0' b <Z>' c + V a , 



whereas at first order, we obtain 

d 2 7 a a + e- 2A Ua a + dA'd z a a + 2g a hr & b d z a c + 



r)V a 

d d g a bc <£' b <Z>' c a d - — a c 



& a (c + d z b) -2V a b = 0. 



Defining a "background covariant" derivative as 

D t <p a = d z + g° 6c <&V, 

and using (fHTj) . we can write ([82]) as 



D 2 + dA'D, + e 



-2A 



□ 



a a -(Vfc-n a bcd $ ,b & d )a c 
& a {c + d z b)-2V a b = 0. 



(81) 



(82) 



(83) 



(84) 
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C.3 Expanding Einstein's Equations 

Einstein's equations are 



= -iV + 2G ab (d^ a )(d^ b ) + -j^9p,V = 0. (85) 

In terms of four-dimensional quantities, this can be written as (normal, 
mixed, and tangential components) 



(n/q)(n/^) - {nicy + n 2 R - An 2 V + 2G ab [(c\$ a )(c>^ b )- 
2n i (d i $ a ){d z $ b ) + (nV - nV)(d;$ a )(d,$ 6 ) 



di{nK?A - V,(n/Cf ) - nK?& Inn + nKldj In 



n- 



(86) 



(87) 



2G ab {d z <S> a -rjdj$ a )d& b = 0, 

and 

-0*(n/C§) + n k V k (nlCi) + n/q(n/C£ + d z \nn- n k d k lnn) + 
nV^j-n + rc/qVj-n* - nlC]V k n i - n 2 R)+ 

An 2 V 

2n 2 G ab (V^ a )(d^ b ) + = 0. 

At zeroth order, the normal and tangential components together 
yield (written on a form suggestive of the superpotential formalism) 

V =l Gab ^-^^lA' 2 , 

2 2 4 (89) 

A" = - —-G ab ^ lb . 
a — 1 

At first order, the normal component gives 

2{d - l)A'c + 4^(L> 2 o a ) - W a a a - 8Vb = 0, (90) 
where we use the notation Q' a = G ab Q' b . The mixed components give 

—not + (d- i)A%b - 2& a d ia a = o. (9i) 

Here, we have used that tj and t)j are transverse, and that Rj = 
— ^e~ 2A de t j and R = at first order. 

d l = to first order (decompose ()9ip as transverse and tangential) , 
and therefore (|91|) implies that 

(92) 



(d - I) A 1 
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Plugging into (|90|) gives 

8V$'a a 



+ 



2V a a a 



2& a D z a a 



(d-l) 2 A' 2 (d-l)A' (d-l)A r 
Using (f89|) . one can show that 



d z b 



2d<S>' 



+ 



2Vn 



d-1 (d-l)A' (d-l) 2 A' 2 (d-l)A' 
Plugging everything into (|84p finally gives us 

D 2 + dA'D z + e~ 2A U 



(93) 



(94) 



rm ^b^d M^' a V c + V a $') 16V& a <f>' 

V\c - n bed® * + , , ,x + 



(d - I) 2 A' 2 



o c = 0. 



(95) 



D Superpotential from the BPS Equa- 
tions 

For Type A, it is possible to find a superpotential W using the BPS 
equations. We will consider spinors in Type IIB SUSY variations that 
satisfy 



e = ie*, T eip e = T $ -e, T^e = e. 



(96) 



The gravitino variation 5tp x = gives 



-2/-2g-2h+| 

16 



e 29 N f - (e 29 - 4e 2h ) N c 



(97) 



where prime denotes differentiation with respect to p. These equations 
are the same as the ones coming from the dilatino variations with 
4> = 4/. Further, 5?pQ = gives 



h' = i e - 2 (/+ ft ) (>/ 2 Ar c + e 2 (/+ fc ) - S> / 2 N 
while 5ip\ = gives 



~ g > = e -2(/ +S ) / 2(/+fe) _ f 



and, finally, (5^3 = gives 



= — e" 



-2{f+h+g) 



. e 0/2 (4^ _ e 2gj ^ _ 4e 2(/ +/l +fc) + 
,2 (f+h+g) _ e 2(f+k+g) _ e |+2g N ' _ 



(98) 



(99) 



(100) 
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The equation of motion for A gives us an expression for the superpo- 
tential 

w = -\^ = -^T- + 2 a + 2h ' + k ') > ( 101 ) 

where (as above) prime denotes differentiation with respect to p. Using 
the above BPS equations, one arrives at 

W = — e '9-^P+x) r e 6 P +x+| / /_ x + e 2 9 \ Nc + N \ _ 

16 L \ (102) 

4 (1 + e 29 + 2e 9+6p+2x ) 



One can check that this superpotential W reproduces the Type A 5d 
potential V through 

V = \G ab W a W b - ^W 2 . (103) 

Also, the equations of motion derived from the superpotential are 
precisely the BPS equations given above. 



E Numerical Methods 

Suppose that we have a system of n scalar fields satisfying a second 
order linear differential equation, and that the boundary conditions 
in the IR single out p linearly independent solutions, whereas the 
boundary conditions in the UV single out q solutions. A solution 
is completely characterized by evaluating it and and its derivative 
at a chosen point. Therefore, let us form vectors (aiRf{\,d p aiR(i\) 
where different % denote different solutions in the IR, and we have 
suppressed the field index. These are p column vectors of size 2n. 
By evolving them numerically from the IR, we can evaluate them at 
any point we like, and therefore they are functions of p. Similarly, 
we form q column vectors from the UV solutions, (auvu), dp&UVU))- 
The question that we need to answer is whether, for a particular value 
of K 2 , we can find a solution that interpolates between the correct 
IR and UV behaviours. In other words, we want to know whether 
we can find a linear combination of the p solutions in the IR and 
write it in terms of a linear combination of the q solutions in the 
UV. For p = q = n, this is true if and only if the deterimant of the 
matrix formed by putting the IR and UV column vectors next to each 
other is equal to zero. It is convenient to evaluate this matrix at a 
point chosen between the IR and UV. In other words, the linearly 
independent solutions satisfying the boundary conditions in the IR 
and in the UV, respectively, are evolved numerically to a midpoint, 
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where the determinant is evaluated. If it is zero for a particular value of 

K 2 , there is a pole in the correlator. This is the midpoint determinant 

method described in [4"2"] . 

We would now like to generalize this method to include cases where 

p and q are not necessarily equal to n. In such cases, the matrix 

obtained by putting the IR and UV column vectors next to each other 

is not generally a square matrix, and therefore we can not answer the 

question of whether the vectors are linearly independent by evaluating 

a determinant. The method we will use instead is the following. First 

we normalize the vectors (aiRn\, d p aiRu\) and (auv(i), dp&UV(i))- Let 

us denote by X" (z = 1, . . . , p + q, a = 1, . . . , 2n) the matrix formed 

by putting these normalized column vectors next to each other. Then 

we construct an orthonormal basis e 1 a (i = 1, . . . ,p + q, a = 1,..., 2n) 

for the subspace spanned by these vectors. Finally, we project the 

normalized vectors onto the basis and form a matrix Y l , = e l n X a ,. 

j j 

This is now the (p+q) x (p+q) matrix whose determinant we compute 
at a midpoint between the IR and UV. Again, if it is equal to zero, 
there is a pole in the correlator. 
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